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In several experiments the magnitude of the contribution of magnetic impurities to the Kondo re- 
sistivity shows size dependence in mesoscopic samples. It was suggested ten years ago that magnetic 
surface anisotropy can be responsible for the size dependence in cases where there is strong spin- 
orbit interaction in the metallic host. The anisotropy energy has the form AE = iCj(nS) 2 where n 
is the vector perpendicular to the plane surface, S is the spin of the magnetic impurity and Kd > 
is inversely proportional to distance d measured from the surface. It has been realized that in the 
tedious calculation an unjustified approximation was applied for the hybridizations of the host atom 
orbitals with the conduction electrons which depend on the position of the host atoms. Namely, the 
momenta of the electrons were replaced by the Fermi momentum Uf- That is reinvestigated con- 
sidering the fc-dependence which leads to singular energy integrals and in contrary to the previous 
result Kd is oscillating like sin(2fc.FG0 and the distance dependence goes like 1/d 3 in the asymptotic 
region. As the anisotropy is oscillating, for integer spin the ground state is either a singlet or a 
doublet depending on distance d, but in the case of the doublet there is no direct electron induced 
transition between those two states at zero temperature. Furthermore, for half-integer (S > 1/2) 
spin it is always a doublet with direct transition only in half of the cases. 

PACS numbers: 72.15.Qm,73.23-b,71.70.Ej 



I. INTRODUCTION 

There are substantial experimental evidences that the 
amplitude of the Kondo effect due to magnetic impurities 
in metallic samples of limited size are reduced^ but with 
unchanged Kondo temperature. That indicates that not 
all of the impurities contribute in the same way. There 
were early speculations that this reduction appears where 
the sample size is comparable with the Kondo screening 
cloud. This is incorrect as the Kondo coupling is lo- 
cal and the only relevant energy scale to be compared 
with the Kondo temperature is the level spacing of the 
conduction electrons which is e.g. zero for semi-infinite 
samples. Later it was suggested 3 that a magnetic surface 
anisotropy can develop due to the spin-orbit interaction 
in the host metal, which has the form 



H = K d {S z ) 



(1) 



where the constant K,i depends on the distance measured 
from the surface of the sample and S z is the component 
of the impurity spin perpendicular to the surface (see 
Fig. [I}. In those papers^^ it was stated that surpris- 
ingly Kd for large distances is always positive and decays 
with the first power of the distance. That result was 
not questioned in Ref. Recently, one of the authors 
(L.Sz.) has called the attention to an unjustified approx- 
imation in the previous lengthy calculation 3 ^ which can 
be responsible for the very surprising results. That ap- 
proximation was that in the hybridizations of the host 
atom orbitals with the conduction electrons which de- 
pend on the position of the host atoms (see Eq.(3) and 



(9) of Ref. [4|]), the momenta of the electrons were re- 
placed by the Fermi momentum kp. That is even not 
the case in the derivation of the Friedel oscillation 6 . 




FIG. 1: The magnetic impurity at a distance d from the sur- 
face in a metallic host with homogeneously dispersed spin- 
orbit scatterers labeled by n. 



Meanwhile great efforts has been made to derive the 
surface anisotropy by using electronic structure calcu- 
lations. First Szunyogh and Gyorffy calculated the 
anisotropy in semi-infinite Au host for Fe impurities 7 . 
They found that Kd is an oscillating function of the dis- 
tance d and the amplitude falls as 1/d 2 . That was a 
calculation of mean field type and the discrepancy be- 
tween those and the analytical ones was not surprising 
as in the latter the diagrams calculated are beyond the 
mean field approximation. Recently, Szunyogh, Zarand, 
Gallego, Munoz and Gyorffy^ have developed another 
model, where the spin-orbit interaction was placed on 
the d-\eve\ of the impurity instead of the host. They 
considered the Friedel oscillation in the density of states 
nearby the Fermi energy due to the presence of the sur- 
face and the different d-orbitals of the impurity coupled 
differently to these oscillations and that is realized in the 
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oscillating anisotropy decaying as 1/eP. It is interesting 
to note that the Hartree-Fock mean field approximation 
and the diagram beyond that play equal role. We are also 
informed that very elaborate calculations by A. Szilva, L. 
Szunyogh, G. Zarand, and M.C. Muhoa^ are in progress 
where the spin-orbit interaction in the host is considered. 

The relative importance of the spin-orbit interactions 
on the d-level and the host material must be very specific 
for which impurity atom and host metal are considered 
and the final answer can be given only by detailed elec- 
tronic structure calculations. 

The present analytical calculation is focused at the os- 
cillating behavior and the decay rate of the K d func- 
tion. All of the results are obtained in the large dis- 
tance asymptotic region, as the preasymptotic calcula- 
tion would be even more difficult. The consequences of 
the oscillating behavior will be discussed at the end of 
the paper. The main goal is to present correct analytical 
result to be compared in the future with the numerical 
results which may lead to the resolution of the present 
discrepancies. The comparison with the experiments is 
left for the future when the numerical calculation will 
be completed by which the very relevant preasymptotic 
behavior is also achievable. 

The Paper is organized as follows. In Section |TT] we 
present the outline of the problem and call the attention 
to the differences in the calculation compared to the ear- 
lier works^. In Section llTll the integrals with respect to 
the energies are performed which are crucial in obtain- 
ing the correct form of the anisotropy. The consequences 
are analyzed in the Conclusions. Some of the matrix 
elements and further details of the calculations are pre- 
sented in Appendices O and [B] respectively. 

II. THE OUTLINE OF THE PROBLEM 

The magnetic impurity scatters the electron in the 
I = 2 orbital channel and the spin-orbit scattering is also 
restricted to those^. As in Ref. [11] we start with the 
conduction electron propagator leaving and arriving at 
the impurity and in meantime it is scattered by one of 
the heavy host atom due to strong spin-orbit scattering. 
The Green's function has a simple form in the coordi- 
nate system where the impurity is in the origin and the 
scattering atom is on the z-axis at a position R n , which 
is called the local system labeled by n. The Anderson 
modeUS is used for the scattering d-levels of the host 
atom and the spin-orbit scattering is assumed to happen 
on the d-level and that determines the symmetry. Fol- 
lowing Ref. Q the conduction electron Green's function 
in first order of the spin-orbit coupling is 

t~ r k?no , .k'm'(T'\^^n) : 

lU! n - £ k 

W kk ,(R n ) 

i mm 1 

+ 12L j (2) 

4^ iuj n - e k (ioJ n - uj d )(iuj n - Lo d ) iuj n - e k , ' 



where now the fc-dependence in W and the w-dependence 
in the d-level Green function is kept i.e. G^{iu> n ) — 
iui X -u) d wnere w rf — £ d — iA and e d (measured from the 
Fermi level) and A are the energy and width of the d- 
level, respectively and — 2 < m < 2 for the conduction 
electrons. 

In Eq. (|2]) now 

W w {R n ) = 

mm 
era 

XV 2 ( B+(k, k')a- + B~(k, k')a + + B z (k, k')a z \ (3) 

which follows from a similar calculation like in Ref. Q 
and A is the strength of the spin-orbit interaction. 
B ± (k,k') and B z (k,k') are 5x5 matrices in the quan- 
tum number m, having the form 

£w(M') = V(3 + m')(2 - m')v m (k)v m ,(k')S m , m , +1 , 

(4a) 



B mrn>( k > k ') = \/( 3 ~ m'){2 + m')v m {k)v m > {k')5 m ,m> -1 , 

(4b) 

B mm'(M') =™v m (k)v m ,(k')8 m , m ,, (4c) 

where the v m (k) matrix elements given in Appendix lAl 
are the same as in Eq. (13) of Ref. [j. These are combina- 
tions of oscillating functions like sm(kR n ) and cos(fci?„) 
combined with powers like (kR n )~ m ~ n (n = 1, 2, . . . ). 




FIG. 2: The self-energy diagram for the impurity spin. The 
double line represents the spin, the single one the conduction 
electrons. The solid circles stand for the exchange interaction 
and the x labeled by n for the effective spin-orbit interaction 
on the orbital of the host atom at R„. 

The next step of the calculation is the rotation of the 
coordinate system from the n-local one to that one where 
the z-axis is perpendicular to the surface. The angle be- 
tween the z-axis of the old (z n ) and the new (z) coor- 
dinate system is labeled by n . The calculation of the 
spin factor of the self-energy diagram (see Fig. ^ giving 
the anisotropy for the impurity spin is also similar to the 
original one (see Eqs.(21)-(25) of Ref. J4J). 

The average over the positions of the scattering atoms 
Rn and R n i must be performed for the whole volume of 
the sample, separately. For the sake of simplicity the 
continuous limit is applied outside the impurity spin. As 
it was shown in the earlier works^^ in order to get the 
dominant contribution one of n-s is nearby the impurity 
and the other one experiences the existence of the surface 
at large distances. 
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The analytical part of the self-energy diagram Fig. [2] 
now is, however, more complicated as W-s depend also on 
four different electronic momenta and the corresponding 
energies appear in the energy denominators of the elec- 
tron Green's functions. In this way the prefactors also 
depend on the momenta and that plays a crucial role in 
the following. 

For the sake of simplicity we consider the conduction 
electron band with constant density of states po i n the 
energy interval —D<e<D where £ is measured from 
the Fermi energy and we will assume linear dispersion, 
i.e. the corresponding fc-values are k = k F + ^- where 
v F is the Fermi velocity. 

As in a noble metal host like Cu, Ag or Au the ri- 
band is below the Fermi energy, it does not give a new 
singularity in the energy integrals (see Section , thus 
we can replace the d-level propagator by a constant e^ 1 . 

Calculating the contribution of the diagram in Fig. [2] 
we applied the Abrikosov's pseudofermion technique 11 
for the spin. After performing the summation over the 
Matsubara-frequencies we get 



D 



D 



Po J dei J de 2 J de 3 J ck 4 

— D -D -D -D 

1 (1 - n F (ei))n F (e 3 ) 1 
£i - £2 £3 + a; - Ei £3 - £4 

+ (£1 <-> £2) + (£3 <-> £4) + (£1 *-* £2 and £ 3 <-> £4) 



x<^ S Fi(R n ,9 n ,R n ,,9 n r,ki,k 2 ,k 3 ,k4) 



+S 2 F 2 (R n ,9 n , Rn>,6 n/ ;ki, fc 2 , k 3 , fc 4 ) 



(5) 



where u> is the energy of the spin after analytical contin- 
uation and the Fi, F 2 functions given in Appendix [Al are 
defined in the same way as in Ref. 0]. 

As that diagram contains two host atoms, averages 
have to be taken over n and n' . According to our simple 
model^ the anisotropy factor— 4 follows as 

D D D D 

K = ^ J d 3 R n J d 3 R n ,pt J efci J de 2 J de 3 J de 4 

-D -D -D -D 

1 (1 - n F (s 1 ))n F (e 3 ) 1 



£1 - £2 £3+^-^1 £3-24 
+ (£1 <-> £2) + (£3 <-» £4) + (si <-> £2 and £3 «-> £4) 

xF 2 (R n , 6 n , Rn>, n >; ki, k2, k 3 , fc 4 ), (6) 



where a 3 is the size of the volume per host atom. 

Changing the order of the summation over the host 
atoms with the energy integrals, the former can be eval- 
uated in a similar way like in Ref. [4J. 



Eq.(29) of Ref. 4] now reads 

D D D D 

K = Pq J dei J de 2 J de 3 J de 4 

-D -D -D -D 

f 1 (l-n F ( £l ))n f (£ 3 ) 1 
\ £1 - £ 2 £3 + w - £1 £3 - £4 

+ (£1 <-> £2) + (£3 «-> £4) + (£1 <-» £2 and £ 3 «-> £ 4 ; 

00 d 



X 



~a? l J dRn'Rft' [Jl(Rn>Rn'i ki, k 2 , k 3 , £4) 

d r 

+J 1 (R n , Rn>;k 3 , fc 4 , fci, fc 2 )] 



dR n R n / dR n iR n ,J 2 (R n ,R n >;ki,k 2 ,k 3 ,k4 



(7) 



where r is a short distance cutoff in range of the atomic 
radius, and 



7T 

Ji{R n ,R n >]ki,k 2 ,k 3 ,k A ) = (2-7r) 2 J d9 n sm0 n 



x / <ffl n , sin 9 n iF 2 (R n ,R n ,,9 n ,6 n > ;k 1 ,k 2 ,k 3 ,k A ), (8) 


J 2 (R n ,R n >;ki,k 2l k 3 ,k 4 ) = (2tt) 2 J d9 n sin9 n 

8n,min 

7T 

x / d9 n > sin F 2 (i? n , i? n ' , 9 n , # n ' ;fci,fc2,fe3,fe4), (9) 



where n , m i„ = arccos(d/ R n ), 9 n ',min = arccos(d/ R n >) 

Since according to the earlier works^ 3 -^ the largest con- 
tribution comes from the first part of Eq. (|7|) correspond- 
ing to Ji we will consider that. 

The evaluation of the integrals with respect to 9 n and 
9 n > gives 

j / r> r> 7 7 / J \ 4 j2 (2KV 2 \, 2 d{R 2 n -d 2 ) 

Ji(R n ,R n ,;k l7 k 2 ,k 3 ,k 4 ) = —J 2 { — ) 3 1 

15 £ it„ 

x [3v (k 2 , R n )vi(ki, R n ) + 3« (fci, R n )vi(k 2 , R n ) 

-2vi(ki, R n )vi(k 2 , R n ) + 2v 1 (k 2 , R n )v 2 (ki, R n ) 

+2wi(fci, R n )v 2 (k 2 , R n ) - 8v 2 (ki, R n )v 2 (k 2 , R n )] 

X [3u (fc4,-Rn') w l( fc 3,-R«') + 3w (fc3,-Rn') w l( fc 4,i?n') 
+Ul(fc3,i?„/)wi(fc 4 ,i?„/) + 2wi(fc 4 ,_R„/)u 2 (fc3,^n') 

+2u 1 (fc 3 , R n ,)v 2 (k 4 , R n ,) + 4v 2 (k 3 , R n ,)v 2 (k 4 , R n ,)] . (10) 



4 



If k%, k%, k 3 , &4 are replaced by k F that gives back the 
half of Eq. (B2) of Ref. [|. 

After a straightforward calculation of the integrals 
with respect to R n and R n > (see Appendix [B]) the first 
part of Eq. (JT]) corresponding to J\ reads 

D D D D 

Pq J dei J de 2 J de 3 J de A 

-D -D -D -D 

1 (1 - n F {e 1 ))n F (e 3 ) 1 
£i - £2 £3 + CJ - £1 £3 - £4 

+ (£1 <-> £2) + (£3 <-» £4) + (£1 <-> £2 and £3 «-> £4) 

4 ,,2ttV 2 A-2 
X - J (— — ) 



15 V £ 







x [C(feid, k 2 d)D{k 3 d, k 4 d) + C{k 3 d, k 4 d)D(k!d, k 2 d)] , (11) 



where the functions C and D are given by Eqs. (|B5[) and 
(jBlip . respectively. 

As C and D are symmetric in their variables we can 
change the integration variables according to the changes 
indicated in the energy dependent factor in the integrand 
of Eq. (|11|) resulting in a simpler form like 

d d D 



Pq J ds\ J de 2 J de 3 J de 4 



-D -D -D 



1 1.4 ^,2irV 2 \, 



4—J z 



£1 - £2 £3 + ^ - £1 £3 - £4 J 15" V £q 
x [C(hd, k 2 d)D(k 3 d, kid) + C(k 3 d, fc 4 d)D(fcid, k 2 d)] ,(12) 

where we have exploited the 1— n F (ei) and ^(£3) factors 
as well. 



III. THE ENERGY INTEGRALS 

In the following the asymptotic behavior for large dis- 
tances d is considered, therefore, only the leading order 
in ^ is kept everywhere. 

For large distances the radial electronic wave functions 
are fast oscillating as the energy is changed. These fast 
oscillations lead to essential cancellations. In order to 
keep track of the cancellations in the limit d — > 00, the 
Riemann theorem with the first asymptotic correction is 
applied in the following fornJ^ 



, „, , . f(b)sm(xb) f(a)sm(xa) 
dsf(s) cos(xs) - — — - — ^ — - (13a) 



a 

and 



' f(a)cos(xa) f(b)cos(xb) 
dsf(s)sm(xs) ~ (13b) 



which is valid in the leading order in 1/x where / must 
be integrable. 

Let us consider the integrations with respect to the 
energies. In the first part of Eq. (TT2"|) the integral with 
respect to e 2 is 



D 



k F + - 



Cjkid, k 2 d) 
£1 - £2 



dk C{k ld ,k 2 d) 
k\ — k 2 



-D 



Introducing a new integration variable Ak = k 2 — k\ 
and using linear dispersion e = v F {k — k F ), the integral 
reads 



d(Ak) 



Cjkid, (fa + Ak)d) 
Ak 



rf 3 



d(Ak) 



x [/ c+ (Jfeid, (fa + A/c)d)(cos[2fcid] 



cos[Afcd] 
Ak 



. r„, , n sin[A/cdl. 
■sin[2M . , J ) 



Afc 



+f c -(hd, (fa + Ak)d) 



cos[Akd] 
Ak 



+f s +(hd, (h + Afc)cO(sin[2M] 



cos[Akd] 
Ak 



r„, ,,sin[Afcdl. 
■cos[2fcid] [ Af J ) 



Ak 



+df s -(k 1 d, (fa + Ak)d) sin[Afcd] 
+Akd 2 f ci (k 1 d, (fci+Afe)d) 

x (Ci[(2ife 1 + Ak)d] - Ci[Akd})] , 



(15) 



where we used the symmetry property of the cosine 
integral function Ci[— x] = Ci[x] and trigonometrical 
identities^. 

To evaluate the terms containing the Ci[x] function we 



use 



Ci[x] 



du- 



dv C ^ (16) 



and change the order of the integrations with respect to 
v and Ak. 

Due to the cosine and sine functions in the integrand, 
the integral is determined by the singularity at Afc = 
(k 2 = fci). Searching for that we expand the /-functions 
around kid in their second variables and then drop the 
terms which are not singular at Afc = 0. Then the inte- 
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gral Eq. jjTSf is 

-kl+kp + - 



a 



d 3 r 

- J d(Ak) 

-fci+fep-JS. 

cos[Afcd 



(/ c +(fcid, fcid) cos[2fcid] 

+/ c _(fcid, fcid) + f s+ (kid, kid) sin[2fcid]) 
sin[Afcd] 



Afc 



-f c+ (kid, kid) sin[2/cid] 



-f s +(kid, kid) cos[2fci<f]) 



(17) 



The range of the integrations can be extended to 
— oo — ► oo as those integrals are independent of d, while 
the added parts are fast oscillating and, therefore, they 
are as it can be proved by using the Riemann 

theorem given by Eq. (fT3|) also. Then using 



d(Ak) 



cos[Akd] 
Afc 



and 



„ . , , sinfAfcdl 



(18) 



(19) 



we get for Eq. ([T 

D 



dei 



Cjkid, k 2 d) 



— D 



d 3 

d 3 



f c +(kid, kid) sm[2kid] — f s +(kid, kid) cos[2fcid] 



3825 
4(fcid)' 
d 3 225tt 
"a3 2(fcid) 5 



sin[2feid] 



cos[2fciG?], 



225 
2(fcid)'' 



■ cos[2kid] 



(20) 



where we kept only the leading order contribution in 
1/kid as kid > k F d 3> 1 according to the range of the 
integration with respect to ki (ei) in Eq. (fT2"j). 

Let's turn to the integration with respect to Si in the 
first part of Eq. (fT2"|) i.e. to 



dei 



£3 + UJ - El 
-D 



ds 



d 3 



dei 



a" j e 3 
o 



C{kid, k 2 d) 

Sl-El 



225^ 1 „ 
-- cos[2fcidJ 



rf 3 



■LO-Ei (kid) 5 2 
dki 



h + £ - *i (hd) 5 2 



2257T 1 r „, „ 
— cos[2fcidJ. 



(21) 



As e 3 < (fc 3 < k F ) in Eq. (JT2J) and u> w the in- 
tegrand has no singularities in the range of the integra- 
tion, thus in order to find the leading order contribution 
in \/k F d we can apply the Riemann theorem given by 
Eq. |[TS]). Then Eq. flS!) is 



dei 



e 3 + uj — ei 
o -D 



ds 



C(hd, k 2 d) 



£l - £2 

d 3 v F 225tt lsin[2fc F d] 
a 3 e 3 + ui (k F df 2 2d 

225?r e F lsin[2fc F d] 
~ (k F a) 3 e 3 + uiA (k F d) 3 



(22) 



in leading order in l/k F d, where we kept only the con- 
tribution of the lower limit, as the contribution of the 
upper limit of the integral in Eq. (|22[) is proportional to 
1/D, thus it is of lower order. 

Now we have to evaluate the remaining integrals with 
respect to £3 and £4 in the first part of Eq. (|12jl i.e. 



\e 3 ^- L 4 M, (23 ) 

£3 + UJ J £3 - £4 

D -D 



Starting with the first part of Eq. (f2"3")) corresponding to 
the g functions in D (see Eq. (IBlip ). we introduce again 
new integration variable Afc = fc4 — fc3 and use linear 
dispersion. Thus Eq. (|2"3"|) reads 



dk 3 



h-k F + ^ 



dk. 



D(k 3 d, kid) 
k 3 - fc 4 



dk 3 



k 3 ~k F + ^ 



-k 3 +k F + - 



d{Ak) D(k 3 d,(k 3 + Ak)d) (M) 



-fe 3 +fc F --2- 



and in the first part containing the g functions (see 
Eq. (jBlip ) we can repeat the considerations used in per- 
forming the integrals with respect to £1 and e 2 giving 



d 3 



dk 3 



k 3 - fcf 



g c+ (k 3 d, k 3 d) sin[2fc 3 d] - g s+ (k 3 d, k 3 d) cos[2fc 3 d] 



•(25) 



Since k F — > we can again apply the Riemann 
theorem given by Eq. (|13[) . The terms coming from the 
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lower limit of the integrals are less by 1/D than the terms 
coming from the upper limit of the integral which give 



d 3 



v F 



-— 7T- 

a 6 u> 



, , ,.cos[2fc F d] , , ,.sin[2fc F <2] 

g c+ (k F d,k F d) — h g B+ {k F d,k F d) 



2d 



2<7 



d 3 v F 



225 cos[2fc F d] 1125 sin[2fc F d] 



2(k F d) 4 2d 
1 e F 225?r 
(k F a) 3 uj 4(k F d)' 



{k F df 
cos[2A; F <i], 



2d 



(26) 



where only the leading order contribution in l/kpd was 
kept. 

Turning to the second part of Eq. (f2"3")l corresponding 
to the h functions in D (see Eq. (|B11I) ). after using the 
following properties 



1 



h VT-,T-,k F r ) 



(27) 



for the h 



c+/c-/s+ 



functions (see Eq. (|B13|) ) and 



(k 3 - k 4 )dh s -(k 3 d, k 4 d, -j) 



kw k F (k F d) 3 



1 h ( k3 kA h \ 
k F k F 



F K F 

for the h s - function, and introducing the s = 
t = ^ new integration variables, we get 

1 Pi(fc F r ,w) 



(k F a) 



15tt 



(28) 
and 

(29) 



where 



P 1 (^,W) = 157T J 



1 + 7 



1 + ^- 



H(s, t, x) 



and 



(30) 



H(s,t,x) = h c+ (s, t, x) cos[(s + t)x] 

+ h c - (s, i, x) cos[(s — t)x] 

+ h s +(s, t, x) sin[(s + t)x] 

+ (s — t)h s -(s, t, x) sin[(s — t)x}. (31) 

Thus the terms corresponding to the h functions give d- 
independent contribution, therefore Eq. (|23p is Eq. (12T)1) 
in leading order in l/k F d. 

Combining that with Eqs. ((22)) and (p~2|) we get for the 
first part of Eq. (fT2"|) in leading order in l/k F d 

. o A 2 A 2 1 ,, . sin[2fc F dl . . 

16e F (Jpo) 2 ^^-, ^Pi(fc F r ,^)- L / J , (32) 



4 {k F af 



(k F df 



where A = irV 2 p$ is the width of the d-levels due to 
hybridization^. 

Turning to the second part of Eq. (fT2"|) , after changing 
the integration variables as s\ <-» £3 and £2 <-> £4 and 
performing similar calculation as before we get in leading 
order with respect to l/k F d 



2 \2 



16e F (Jp ) 



where 



,A 2 A 



£q {k F a)' 



sin[2fc F d] 

■^^'^'IfcwF' (33) 



1+7 



1+7 



P 2 (x,u) = 15tt J 



ds- 



1 - 



dt 



H(s, t, x) 



Thus the anisotropy factor coming from the first part 
of Eq. (JTJ) corresponding to J% - giving the leading 
contribution^ 3 -^ - in leading order in l/k F d is 



2 \2 



K = 16e F (Jp ) 



,A 2 A 



1 



£0 (fc F a)' 



: P(/c F r , u) 



sm[2k F d] 
(k F d) 3 



where 



P(x,u>) = Px(x, co)+P 2 (x,uj) 



(35) 
(36) 



which is finite for w = 0. 

Evaluating P by numerical integration it turns out that 
for x ~ 1— the integrals with respect to t and s are dom- 
inated by the t = s and s — 1 singularity, respectively. 
Thus 



P(x, W = 0)«157r J d S ^^ x) 



In 



s - 2 



15irH(l,l,x) 



(37) 



where we used D — e F . In Table U we compare the results 
obtained by numerical integration and by Eq. (|37p for 
a; ~ 0.5 -1.5. 



X 


0.5 


0.75 


1 


1.25 


1.5 


P num (i,w = 0.01) 


-135.8 


-428.1 


-916.5 


-1552.9 


-2215.9 


P appr (x,w = 0) 


-138.3 


-438.4 


-952.6 


-1664.6 


-2514.1 



TABLE I: Comparison of P num (a;,cj = 0.01) obtained by nu- 
merical integration and P appr (a;,tj = 0) obtained by Eq. (|37p 
for x ~ 0.5 — 1.5. 



The result Eq. (f35|) for the anisotropy factor is es- 
sentially different from the earlier one (see Eq. (31) of 
Ref. Q) obtained by an unjustified assumption which 
corresponds formally to the approximation C{k\d, k 2 d) « 
C(k F d, k F d) and D(k 3 d, k 4 d) w D(k F d, k F d) in Eq. (fT2|) . 
The main differences are 

• it contains the oscillating factor sin[2fc F d] 
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• the asymptotic distance dependence is rer^gp in- 
stead of jt^j, thus it is essentially weaker 

• instead of the /(^)P oW (fc_F?'o) factor which con- 
tains the short range cutoff r$ and estimated to be 
between 50 — 950 for w^O and fc^o ~ 0-5 — 1.5, 
we have P new (kpro,uj) which is for u> = and 
k F r ~ 0.5 - 1.5 between -140 and -2500, 

thus in the asymptotic region using the same parameters 
as in Eq. (32) in Ref. [3 



0.01 



cV < \K d \ < 



1.75 



■eV. 



(38) 



Finally, the question can be raised how justified is the 
assumption of homogeneous distribution of the spin-orbit 
scatterers. In order to give the answer in the following 
the scatterers are considered homogeneously distributed 
on sheets which are parallel to the surface and separated 
by a distance a. According to the previous works 3 ^ the 
pair of sheets in equal distances from the impurity do 
not contribute to the anisotropy, thus only the unpaired 
ones must be considered. The sheet n is in the distance 
na from the impurity and only sheets with n > d/a are 
considered. As it was discussed in Ref. 

ma 

one of the 

two scatterers n and n' is nearby the impurity and the 
other one is far from it on one of the sheets considered. 
Therefore, the contribution of the sheets are additive. 
The contribution of sheet n can be easily obtained from 
the present calculation as 



K„ 



9 K 
a dd K ^ 



(39) 



where the derivative gives the contribution of an infinitely 
narrow layer and the prefactor provides the correct nor- 
malization. Thus the final result in the asymptotic region 
is 



K„ 



E 



2ak F 



cos(2fcpna) 
(k F na) 3 



(40) 



where the omitted prefactor is the same as in Eq. (|35[) . 

Thus the separate sheets contribute by different signs 
and amplitude. Due to the fast decay by increasing n 
only sheets of restricted numbers are essential. Adding 
the contribution of the sheets with different signs and 
amplitudes (likely randomly distributed) the final ampli- 
tude of the anisotropy \K | can be larger than \K d \ thus 
the decay rate can be somewhat reduced. The situa- 
tion is different in a coherent case with kpa = pn where 
p is integer. For even p the contributions have the same 
sign and ^ ^ ~ §^ which provides a slower decay rate. 



IV. CONCLUSIONS 

The amplitude of the anisotropy is oscillating and 
weaker than the one earlier estimated 3,4 . In these 



changes the sharp edge in the i?-integral thus the ex- 
istence of the surface is crucial. We used a uniform dis- 
tribution of the spin-orbit scatterers in space. Consider- 
ing the question how the results are changed in the case 
where continuum layers of the scatterers are considered, 
it is argued that the overall behavior is not expected to 
change, but only the amplitude can be influenced mod- 
erately. Furthermore, it is assumed, that the spin-orbit 
scattering is point like, but a finite extension r d (it is 
assumed that k F r d < ir) smears somewhat the sin[2/sprf] 
function in Eq. (|3"5"|) . The actual size of that can be esti- 
mated only by electronic structure calculations and cer- 
tainly will be included in the recent work of one of the 
authors (L.Sz.) and his coworkers in progress 9 . 

The way how the spin is frozen by the surface 
anisotropy is essentially different from the previous 
works^ as K d is not always positive (see Fig. [3]) and 
also depends whether the spin is integer or half-integer. 



K d >0 



K d <0 



S z =±2 
S z =±l 
Sz-0 



singlet 



K d >0 



S z =±5/2 

S z =±3/2 
S z =+l/2 



S z =0 
S 7 =± 1 



no scattering 
(a) 



S z =±l/2 
S z =±3/2 
S z = ± 5/2 



doublet 



K d <0 



doublet 
scattering 



doublet 
no scattering 



(b) 



FIG. 3: The level splitting due to the surface anisotropy for 
(a) integer (e.g. 5 = 2) and (b) half-integer (e.g. 5 = 5/2) 
spins. It is also indicated whether the electrons can be scat- 
tered by the degenerate ground states or not. 



• Integer spin (e.g. S = 2): the ground state is ei- 
ther singlet or doublet depending on the sign of 
K dl however, the electron cannot be scattered by 
transition in the doublet as the spin momentum 
difference is AS = ±4 and the turning the electron 
spin allows only AS = ±1. Thus the spin at low 
temperature can be completely frozen in. 

• Half-integer spin (e.g. S — 5/2): the ground state 
is always doublet, however, in one of the cases 
AS = ±1 which can cause scattering in contrary 
to the case where AS = ±5. Thus only half of 
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the impurities can cause electron spin flips at low 
temperature. That is different from the previously 
assumed K~d > case. 

This result shows analogies with Ref. [1 51 ] where mag- 
netic molecules with large spins on metallic surface were 
considered. The spin levels are split in a similar way and 
electron induced transitions are allowed only between cer- 
tain levels. 

For comparison with experiment the preasymptotic be- 
havior is very essential which is beyond the scope of the 
present paper. The electronic calculations in progress^ 
must provide those information including the amplitude 
of the anisotropy as in the other model in Ref. [8j . 

The present results valid in the asymptotic region can 
provide a good possibility to test the numerical calcula- 
tion. Detailed comparison with experiments must wait 
for completing the numerical calculation which is going 
to provide more necessary information. 



where 

&2 &s (74 V W klk2 (R n ,e n )W kik3 (A3) 

' J m\rri2 m 2 m 1 

mm' <7± &2 £74 0"3 

W gives the form of W in the rotated coordinate system 
(where the z-axis is perpendicular to the surface) given 
by 

W kk' (RniOn) = &rn+<7,m'+a' ^ d^^{d n )d^ 2 \d n ) 

mm _ _ / 

/ mm 

00 / 

aa 

■w ky (Rn)d$ ml (-e n )d%*\_ en) (A4 ) 

mm' 

where the Wigner-formula for rotation matrices^ was 
used. 
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APPENDIX A 

The v m (k) matrix elements are the same as in Eq. (13) 
of Ref. 1: 



v (k,R n ) = 10 



sin(fci?„) 3cos(fci?„) 12sin(fci?„ 



2kR n (kR n ) 2 
27cos(fci?„) 27sin(fci? n ) 



(kR n y 



{kR n f 



(kR n ) 3 
(Ala) 



vi(k,R n 



15 



(kRn) 2 
12sin(fci?„) 
(kR n ) 5 



(kRn) 3 



(kRn) 4 
(Alb) 



v 2 (k, Rn) = 15 



sin(fci?„) 3cos(fci?„) 3sin(fci?„) 



(kR n f (kRn) 4 (kR n ) 5 

(Ale) 

The F\ and F 2 functions are defined in the same way 
as in Ref. [4J, namely 

2J 2 

F 1 (Rn,9n,Rn>,0n>;ki,k 2 ,k 3 ,k4) = —7-/1122, (A2a) 



and 



F 2 (R n , On, Rn',0 n ';ki, k 2 , ks, ki) 



2 J 2 

— (/mi - /1212 - /1122), (A2b) 



Here we perform the integrations with respect to R n 
and R n ' in the first part of Eq. (fTTj) . Let us start with 
the integration with respect to R n , namely 



C(ki,k 2 ) := -1 / dR n R 



,2 d(Rl - d 2 ) 



R 3 

x [3v (k 2 , R n )vi(k u R n ) + 3v (ki, Rn)vi(k 2 , Rn) 
-2vi(ki, R n )vx(k 2l R n ) + 2vi(k 2l R n )v 2 (ki, R n ) 
+2v 1 (k 1 ,R n )v 2 (k 2 ,R n )-8v 2 (k 1 ,R n )v 2 (k 2 ,R n )](Bl) 

After substituting the v m (k, Rn) matrix elements given 
by Eq. (|Aip and introducing the dimensionless integra- 
tion variable y = R n /d and notations t\ = k\d and 
t 2 = k 2 d we get 



C(k 1 ,k 2 ) = C(k 1 d,k 2 d), 



(B2) 



cos(kR n ) 5sin(fci?„) 12cos(/ci?„) where 



d 3 

C(h,t 2 ) = \ 



dy- 



h)y] 



x [C c+ (ti,t 2 ,y) cos[(h 
+C c -(ti,t 2 ,y) cos[(*i - t 2 )y] 
+C s+ (t 1 ,t 2 ,y) sin[(ti +t 2 )y] 
+C s -(t 1 ,t 2 ,y) sin[(ti -t 2 )y]]. 

The occurring integrals with respect to y look like 

00 

y 2 -1 cos[ty] 



and 



G c (t,n) 



G s (t,n) 



dy- 



y 



y h 



(B3) 



(B4a) 



dy 



y 2 - 1 sm[ty] 



(B4b) 
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Evaluating the G s and G c functions analytically by 
using MATHEMATICA, we get 

d 3 

C(t u t 2 ) = -z[f t *(t U t 2 )cOB[t l + t2] 

a 6 

+fc-(h,tz) cos[ti - t 2 ] 
+f s+ {t 1 ,t 2 )sm[t 1 +t 2 ] 

-*a)/—(*i,*a)siii[ti -*a] 
+(ti-t2) 2 fci(ti,t2)(Ci[t 1 + t2}-Ci[t 1 -t 2 })],(B5) 

where for t\ = t 2 = t 3> 1 

3825 



fc+(t,t) 
fc-(t,t) 
fs+(t,t) 

/.-(M) 
fci(t,t) P 



4i 6 
225 

^ 225 

^ 225 
* 2F 
225 

2F 



(B6) 



and Ci[t] — — J du^ 2 ^ is the cosine integral function. 
t 

Let us consider now the integration with respect to R n 
in the first part of Eq. (TTTj) : 



£>(h,k 4 ) := i / <//?„-/?:, 



x [3t;o(fc 4 , R n >)vi(k 3 , R n >) + 3v Q (k 3 , R n >)vi(k 4 , R n >) 

+V 4 {k 3 , R n ')vi(k 4 , R n ') + 2vi(k 4 , R n ')v 2 (k 3 , Rn') 

+ 2 Vl (ks , Rn' )v 2 (h,R n ')+ 4w 2 (fc 3 , i?„' )«a (^4, -Rn' )] (B7) 

After substituting the v m (k, R n ') matrix elements given 
by Eq. (|A1|) and introducing the dimensionless inte- 
gration variable y' = R n >/d and notations £3 = k 3 d, 
t 4 = k 4 d, and yo = r o/d we get 

L>(fc 3 , fc 4 ) = D(k 3 d, k A d), (B8) 

where 

1 

I>(*3,*4) = ^3 / W 2 



x [D c+ (t 3 , t 4 ,y') cos[(i 3 + t 4 )y'] 

+D c _{t 3 ,U,y')cOs[{t 3 -U)y'] 

+D s+ (t 3 ,t 4 ,y')sm[(t 3 +t 4 )y'} 
+£>._(t 3 , *4, y') sm[{t 3 - t 4 )y']] ■ (B9) 
The occurring integrals with respect to y' look like 

1 

~ rrvzl+n \ 

(BlOa) 

J/o 



. , n cos[V] 
H c {t,n) := / y — . 



and 



H s {t,n) 



dy'y-^f 



(BlOb) 



Evaluating the H s and H c functions analytically by 
using MATHEMATICA, we get 

d 3 

D(t 3 ,t 4 ) = — \g c+ (t 3 ,t 4 ) cos[t 3 + t 4 ] 

+g c -{t 3 ,t 4 )cos[t 3 - t 4 ] 

+# s+ (t 3 ,t 4 )sin[i 3 + i 4 ] 

+(*3 - ti)9s- C*3> sin[£ 3 - t 4 ] 

+/i c+ (i 3 , i 4 , ?/o) cos[(i 3 + t 4 )j/ ] 

+h c _(t 3 ,t 4 ,y ) cos[(i 3 -t 4 )j/ ] 

+/i s+ (t 3 , i 4 , j/o) sin[(t 3 + t 4 )yo] 

+(t 3 - h)h s - (t 3 , t 4 , y ) sin[(t 3 - * 4 )t/ ]] , (Bll) 

where for t 3 — t 4 = t 3> 1 



and 



h c +(t 3 ,t 4 ,y ) 



h c -(t 3 ,t 4 ,y ) 



h s+ (t 3 ,t 4 ,y ) 



and 



h s ~{h,t 4 ,Vo) 



9c+{t,t) 


225 




9c-{t,t) 


^ 225 
~ ~2t* 




9s+{t,t) 


1125 

t 5 




9s-(t,t) 


1125 

2t 6 




20250 


8100 


20250 


tot a Un 


Lob a Un 
4 t?U 


tot a Ur\ 


8100 


1350 


6525 


4£*U 


*3*42/o 


2*3*42/0 


1350 


225 




+ +4+2 3 

t 3 t 4yo 


2*|*lz/o ' 




20250 


8100 


20250 


J 3*42/o 


+ tMvo 


*3*42/o 


8100 


1350 


6525 


+ my§ + 


*i*42/o 


2*|*42/o 


1350 


225 




*t*42/o 


2*§*l?/o ' 




20250 


20250 


1350 


*t*42/o 


*3*42/o 


*3*42/o 


8100 


8100 


1350 


*3*42/o 


* 3 *4 vt 


*3*42/o 


1125 


1125 




2*3*42/0 


2*3*42/0 


20250 


1350 


6750 


*l*42/o 


*3*42/o 


*3*4?/0 


1350 


1125 




*i*42/o 


2*3*42/o 





(B12) 



(B13a) 



(B13b) 



(B13c) 



(B13d) 



!/(> 
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